Reduced phase error through optimized control of a superconducting qubit 
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Minimizing phase and other errors in experimental quantum gates allows higher fidelity quantum 
processing. To quantify and correct for phase errors in particular, we have developed a new experi- 
mental metrology — amplified phase error (APE) pulses — that amplifies and helps identify phase 
errors in general multi-level qubit architectures. In order to correct for both phase and amplitude 
errors specific to virtual transitions and leakage outside of the qubit manifold, we implement "half 
derivative" an experimental simplification of derivative reduction by adiabatic gate (DRAG) control 
theory. The phase errors are lowered by about a factor of five using this method to ~ 1.6° per gate, 
and can be tuned to zero. Leakage outside the qubit manifold, to the qubit |2) state, is also reduced 
to - 10" 4 for 20% faster gates. 



Many candidate systems for quantum computation dis- 
play several quantum energy levels, with computing ar- 
chitectures either employing just two of these levels [1- 
12], or using three or more in ("d-level") qudit-b&se& ap- 
proaches [2, 6, 13]. Achieving good control over the com- 
putational states, while avoiding interference from the 
non-computational ones, requires measuring and under- 
standing the source and magnitude of error-generating 
processes as well as correcting for such errors. For ex- 
ample, gate operations in the superconducting phase [2] 
and transmon [6] qubits, which are relatively weakly an- 
harmonic, can generate leakage from the computational 
manifold to higher (non-computational) levels ("mani- 
fold leakage"). However, by careful experimental design 
one can identify and minimize the source of such errors, 
thereby improving qubit control. 

Previous experiments identified the amplitude errors 
associated with manifold leakage by using error budget- 
ing and metrology [14]. By directly measuring leakage 
using the Ramsey error filter protocol, this error is now 
understood and can be suppressed to ~ 10 -4 [14], at 
or near the threshold for fault-tolerant quantum com- 
puting [15, 16]. Randomized benchmarking, an alter- 
native approach to measuring gate error, relies on a 
many-pulse protocol, effectively averaging over the Bloch 
sphere to quantify gate fidelity [16, 17]. Unfortunately, 
this approach optimizes a single value (gate fidelity) and 
does not distinguish between amplitude and phase er- 
rors. Quantum process tomography by contrast provides 
a complete analysis of gate operation [1, 18, 19], but re- 
quires calibrated X^^ and Y^/i pulses, which can them- 
selves be error sources. 

Phase errors also contribute to overall gate error, but 
the relative contribution differs from amplitude-related 
errors. Consequently, a method is needed to separately 
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quantify the phase error generated by a gate. This will 
aid in identifying the source of these errors, improve cal- 
ibration of high-fidelity tomography pulses, and provide 
a benchmark for optimizing control pulses. 

In this Letter we explicitly show that gate errors can be 
separated into amplitude and phase errors. Virtual ex- 
citations to non-computational levels (e.g. the |2) state) 
during gate operation result in phase errors, while real 
excitations after the gate produce amplitude errors. In 
order to better quantify the phase errors, we introduce 
a new metrology, "amplified phase error" (APE) pulses, 
which uses a Ramsey fringe experiment to measure and 
amplify this ubiquitous source of error. We focus on er- 
rors related to 7r/2 pulses, because such pulses provide 
the basis for tomography and are essential in algorithms. 
We also demonstrate a simplified experimental version of 
the protocol [20] termed "derivative removal by adiabatic 
gates" (DRAG) [15], which we call "half-derivative" (HD) 
pulses. By using HD pulses together with APE metrol- 
ogy, we measure and reduce the phase error to 1.6° per 
gate, a factor of five reduction from unoptimized perfor- 
mance. As a demonstration of this method, we perform 
quantum state tomography to map out the trajectories 
of typical HD pulses, including a 7r-pulse and a rotation 
about an arbitrary axis, using X, Y and Z controls to im- 
plement an (off-equator) Hadamard gate. In addition, we 
show that APE metrology is a universal tool for probing 
phase errors on any of the higher qudit levels. 

In the experiments described here, we used a sin- 
gle superconducting phase qubit with 7\ = 450 ns and 
T| cho = 390 ns. The circuit layout and operation have 
been described previously [2]. We have three-axis control 
over the qubit: microwave pulses of arbitrary amplitude 
and phase, resonant with the qubit |1) o |0) transition 
frequency /io, produce rotations about any axis in the 
x — y plane, while current pulses on the qubit bias line 
adiabatically change the qubit frequency, causing phase 
accumulation between |0) and |1), generating z-axis rota- 
tions [21]. The |2) o |1) transition frequency /21 differs 



from /io by A/2tt = / 2 i -/io ~ -200 MHz (see Fig. la). 

Non-ideal qubit behavior can arise from both leakage 
at the end of the gate and virtual transitions to higher 
states during on-resonant operations. The leakage is an 
amplitude error, representing loss of probability to states 
outside the manifold. Leakage can be reduced to ~ 10 -4 
by careful shaping of the microwave envelope and choos- 
ing the correct gate duration, which scales as 1/A [14, 22]. 

The phase error arising from virtual transitions (espe- 
cially to the 1 2) state) can be modeled as effective qubit 
rotations about the z-axis. We first restrict ourselves to 
simple gates comprising i\ and tt/2 rotations. An X^^ 
pulse (a rotation about the x-axis by an angle 9 = tt/2) 
ideally produces the transformation 
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where cr x is one of the Pauli matrices. From numerical 
simulations of our multi-level qubit, we find that this 
transformation, expressed in quantum circuit language, 
is instead 



X n/2 



Z e X^i2Z e 



(2) 



where Z e is the phase error of interest and < e <C 1 [23] . 
The leading term in Eq. (2) is a global phase and can 
be ignored. We note that Eq. (2) differs from X*, 2 = 
Z_ e X^i 2 Z ei which corresponds to a rotation about a new 
axis e away from the x-axis in the x — y plane. 

The phase error e is a function of both the rotation 
angle and the gate time t g . From simulations, we find 
that e ~ @ 2 /tg- Longer gate times decrease the virtual 
transitions and consequently reduce the phase error, con- 
sistent with the AC-Stark effect [24]. 

In order to best reduce this error, we first sought a 
protocol that would amplify the error e. If we consider a 
27r rotation generated by concatenating four tt/2 pulses, 
we find from Eq. (2) that X'\ 2 = — e~ Ale J, where I is 
the identity [23] . A concatenated 2tt rotation thus does 
not accumulate the relative phase error. 

We next examine the pseudo-identity operation that 
is formed by concatenating positive and negative 6 rota- 
tions. For a first-order expansion with e <C 1 we find 



I' e = (Z £ X Z e )(Z € X_ Z e ) 

'\ + i(cos0 — l)e — (sin#)e 

(sin 6) e 1 — i(cos + 3)e 



, (3) 



where Xq is an arbitrary rotation of 6 about the x- 
axis[23]. For = n we find that i£ = e~ 2ie 7, which 
is similar to the 2tt rotation, as the phase error e cancels. 
However, for = tt/2 we find 



^tt/2 — (Z e X_ 7V /2Z e )(Z e X 7V / 2 Ze) 



1 — ie e 
— e 1 — 3ie I ' 



showing phase error accumulation. 
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FIG. 1: (Color online) Multilevel qubit and amplified phase 
error (APE) metrology, (a) Qubit with three energy levels, 
illustrating phase error due to virtual transitions to the |2) 
state. While performing on-resonance |0) <H> |1) gate oper- 
ations at frequency /io, virtual transitions to |2) create a 
phase change in |1). (b) Left, single-control (X-quadrature) 
APE pulse sequence, where the pulse shape represents the mi- 
crowave envelope. Ramsey-fringe experiments are modified by 
-^Tr/2 ( n ^ {Ojl^S^}) pseudo-identity operations in-between 
the first X^/2 and last 0^/2 pulses. At the end the Z con- 
trol line is pulsed to measure the probability of |1). Right, 
the probability of measuring |1) as a function of the rota- 
tion axis (j) of the final 7r/2 pulse. Fits to extract the phase 
shift are plotted as lines with the data (dots), (c) Left, same 
pulse sequence as in a, but with the addition of y-quadrature 
"half-derivative" pulses, as discussed in the text. Right, the 
data (dots) with fits (lines) show small phase shifts for the 
HD pulse sequence. 



To measure this error, we combine the result from 
Eq. (4) with a Ramsey fringe experiment, forming the 
"amplified phase error" (APE) sequence. The APE se- 
quence consists of inserting n successive 1^ , 2 pseudo- 
identity operations between the tt/2 pulses that define 
a Ramsey fringe measurement (Fig. lb). The phase er- 
ror is amplified by 2n for n applications of the pseudo- 
identity operation, 
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By applying APE pulses to the state \i/j) = (\0)-i\l))/y/2 
followed by a final (p^/2 pulse, we directly probe the phase 
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FIG. 2: (Color online) Amplitude errors due to leakage into 
the 1 2) state from an X^-pulse. Plot of |2) error versus pulse 
width (FWHM) r for single-control Gaussian (black dots) 
and HD (blue squares) 7r-pulses. The 6 ns HD pulse pro- 
duces x5 lower error and 20% faster gates. The solid lines 
are three-state simulations using Gaussian (black) and HD 
(blue) pulses. 



error due to the X^j 2 pulses. 

Figure lb shows the probability of measuring the |1) 
state versus rotation axis <\> of the final 0^/2 pulse, for 
I'™i 2 (n = 0,1,3,5) pseudo-identity operations. Consis- 
tent with Eq. (5), the phase error scales with n. For n = 5 
in Fig. lb, the final pulse is ~ 90° out of phase, corre- 
sponding to a 10 x phase error amplification from a total 
of 11 pulses, yielding 7.3° phase error per gate. The os- 
cillation amplitude is also reduced, due to decoherence. 

To correct the phase error, we employ the derivative 
reduction by adiabatic gates (DRAG) protocol [15]. The 
original DRAG prescription uses three controls, X, Y, 
and Z. The X control provides the original envelope- 
shaping to the microwaves, which we implemented as 
a Gaussian in time with arbitrary amplitude A, X = 
A exp[-41n(2)(t - t ) 2 /r 2 ], where r is the full-width-at- 
half-maximum (FWHM) and to the time at the center of 
the pulse. The quadrature control Y = —X/A is the time 
derivative of the X control scaled by the nonlinearity A. 
The Z control produces a dynamic detuning pulse during 
the gate that removes the effective z-rotations from the 
virtual transitions. 

We find both in simulations and experiment that the 
Y and Z controls are not independent. For experimen- 
tal simplicity, we set Z to zero and compensate by re- 
ducing the magnitude of the Y control by 1/2, to form 
the so-called "half-derivative" (HD) protocol [23]. For a 
Gaussian envelope on the X control, the HD pulses are 
as illustrated in Fig. lc and differ from the DRAG pulses 
by the quadrature controls, Y = — X/(2A), Z = 0. Ex- 
perimentally, we found that HD performed as well as full 
DRAG. 

The HD pulse sequence in Fig. lc is the same as Fig. lb, 
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FIG. 3: (Color online) Gate trajectories using quantum state 
tomography (QST) for single control and HD pulses, (a) 
Left, pulse sequence for an X rotation. A Gaussian pulse 
of fixed length (FWHM = 6 ns) and varying amplitude com- 
pletes the rotation angle to — it. QST is performed at 
every incremental increase of amplitude. Right, the recon- 
structed quantum state data is plotted on the Bloch sphere 
from two perspectives, looking down the x and the —y axes, 
(b) Same as in a except using HD (X and Y simultaneous 
control) pulses to perform the X rotation. No phase error is 
observed, (c) Left, pulse sequence of a two-part trajectory 
using HD pulses and Z control to form a Hadamard gate. X, 
Y, and Z controls are fixed- length (FWHM = 6 ns) pulses 
with increasing amplitudes to execute tt/\/2 rotations about 
X and Z. Trajectory 2 pulses ramp-up only after trajectory 
1 pulses are at full amplitude. Right, the reconstructed QST 
data. Each trajectory completes a Hadamard gate, taking 
|0}^(|0> + |1>W2^|0>. 



with the addition of the Y control. Data are plotted 
for the same number of I\ 2 pseudo-identity operations. 
We find by applying the HD protocol the phase error is 
reduced to 1.6° per gate, and can be further minimized 
by tuning the amplitude of the derivative pulse. 

HD pulses also reduce the leakage to the |2) state. Plot- 
ted in Fig. 2 are the data from a Ramsey error filter [14] 
for both single control Gaussian and HD pulses. A 6 ns 
(FWHM) HD X n pulse gives a |2) state probability of 
10 -4 , almost an order of magnitude better than a non- 
HD pulse of the same width, which consequently provides 
a 20% faster gate [14]. 

With calibrated X n / 2 an d 1^/2 pulses, we can perform 



quantum state tomography (QST) without worry of mis- 
calibrated measurement axes. As a practical demonstra- 
tion of how HD pulses reduce phase error, we perform 
QST [21] with and without HD. Figure 3a(b) shows the 
pulse sequence and data for the Gaussian pulses (HD 
pulses) during an Xq rotation. The pulses are of fixed 
length (FWHM = 6 ns) with variable amplitude 9. QST 
is performed at each incremental increase of amplitude 
and the quantum state is recreated in the Bloch sphere 
as shown to the right of each of the respective pulse se- 
quences. In contrast with the single control Gaussian 
pulses, the HD pulses execute a meridian trajectory with 
no phase error with increasing 0. 

The final HD demonstration is an (off-equator) 
Hadamard gate, shown in Fig. 3c, which uses all three 
control lines [23] . We incrementally increase the am- 
plitude of all three control lines using fixed length 
(FWHM = 6 ns) pulses to perform rotations from to 
ir/\/2 about both the x and z axes, which at full ampli- 
tude gives the Hadamard gate H (|0) -> (|0) + |l))/>/2). 
The trajectory concludes with a second set of pulses to 
complete the identity operation I = HH, and returning 
to the initial state (|0) + |l))/>/2 -> |0). 

We also consider the challenge of optimizing control 
pulses for each d- level of a qudit [2] . Tomography verifies 
the operation, but again relies on the calibrations of 7r/2 
pulses. In principle, APE provides the necessary phase 
calibration certification for tomography. 

To demonstrate the general utility of the APE proto- 
col, we implement this metrology on the qudit level for 
state 1 2). We first calibrate the (7r)io an d (tt)21 pulses 
to generate the |0) — >• |1) and |1) — >• |2) transitions, re- 
spectively [2]. As shown in Fig. 4, we first prepare the |1) 
state via an HD tt pulse so that we can then perform a 
Ramsey fringe experiment using the |1) and |2) states. 
The APE pulses are applied between the first and last 
(71-/2)21 pulses, only now resonant with f 2 \. The data 
are for single controlled Gaussian envelope pulses, i.e. no 
HD protocol for the |1) — )> |2) pulses. 

Surprisingly, after n = 5 pseudo-identity operations, 
only 12° of phase error is measured, equivalent to 1.1° per 
gate. We offer a qualitative interpretation: the relatively 
small anharmonicity of the phase qudit, combined with 
the symmetric virtual transitions to the |3) and |0) states, 
provide complementary phase shifts that partially cancel 
out the phase error. 

In conclusion, we introduced a new metrology tool, 
amplified phase error (APE) pulses, which can amplify 
the phase error by an order of magnitude. Together with 
APE and half-derivative pulses, our simplified variant of 
DRAG [15], we identify and reduce phase errors to 1.6° 
per gate. By simply re-scaling the analytic form for the 
HD pulses, the phase error can be completely removed. 
The HD pulses also can increase gate speed by 20%. 

We acknowledge discussions with Jerry Chow, Jay 
Gambetta, Felix Motzoi, and Frank Wilhelm. Devices 
were made at the UCSB Nanofabrication Facility, part of 
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FIG. 4: (Color online) Multilevel qudit and APE metrology 
for the |1) <H> 1 2) transition. Inset, the four energy levels of 
the qudit illustrate the nearly symmetric virtual transitions 
to states 1 3) and |0) that mostly cancel the phase error during 
the /21 resonant drive. Left, pulse sequence showing (7r)io for 
a |0) — > |1) transition then a Ramsey fringe with the pseudo- 
identity operation on the |1) ++ |2) transition. Right, the 
probability of measuring the |2) state as a function of the 
rotation axis of the final (0^/2)21 pulse for n G {0,1,3,5} 
pseudo-identity operations. With only a single control pulse 
(no HD correction), the phase error is 1.1° per gate. 
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TRACKING QUBIT FREQUENCY 

Optimized qubit control pulses rely on a precise mea- 
surement of the qubit transition frequency. After per- 
forming spectroscopy to find the qubit frequency fio to 
within ~ 1 MHz, we fix the gate time r and tune the 
microwave amplitude to execute a ir rotation. We verify 
that the computed amplitude for tt/2 rotation is indeed 
half that for a tt rotation by performing two consecutive 
7r/2 rotations and comparing the probability P\ between 
that of a single ix rotation. 

Next, to precisely measure the qubit frequency, we use 
a Ramsey fringe experiment, where the final ^^-pulse 
rotates at 50MHz about a variable axis on the equator 
of the Bloch sphere. A frequency shift in the oscillations 
of Pi different from 50MHz is the amount the microwave 
drive is detuned from the qubit frequency. Correcting 
for this offset precisely tunes the microwave drive to the 
qubit frequency to within 1 part in 10 4 (sub MHz reso- 
lution), which is consistent with limits set from 1/f flux- 
noise fluctuations [1]. 

We also can perform a more complete test (2-D scans) 
of this frequency calibration by noting that the phase er- 
ror is e = 5f5t if the microwave drive is 5f off-resonance 
from the qubit frequency for some time 5t. Therefore, we 
verify that the microwave carrier matches the qubit fre- 
quency for the entire St of the APE sequences via a Ram- 
sey fringe experiment as shown in Fig. SI. Fig. Sib shows 
the data for two different detunings: left, Sf = 10 MHz 
and right, 5f < 1 MHz (after performing the calibration 
described above). When the microwave drive is detuned 
by 5 = 10 MHz the data show a distinct tilt and clear 
oscillations with a frequency of 10 MHz. After calibra- 
tion the data has no beating and therefore no sign of a 
detuned microwave drive. This confirms that the qubit 
frequency is tracked precisely throughout the duration of 
the desired (Amplified Phase Error) sequences. 

To compute the nonlinearity A/2tt = /21 — /10 f° r use 
in the HD protocol, we directly measure the transition 
frequency between states |1) and |2) with a Ramsey Error 
Filter (REF)[2]. The REF uses an oscillation, provides 
finer resolution, and is simpler to automate experimen- 




0.24 



[rad] 



FIG. SI. (Color online) Tracking the qubit frequency fio. (a) 
Ramsey fringe sequence that consists of two Gaussian-shaped 
7r/2-pulses, separated in time by t, followed by a measure 
pulse on the Z-control tuned to tunnel the |1) state. The first 
7r/2-pulse defines the rotation axis; by convention this is the 
x-axis. The second pulse is delayed by a time t with variable 
rotation axis (f). (b) Plots of Pi versus separation time t and 
phase 4> for pulse sequence in a. Left, microwave drive is 
10 MHz detuned from the qubit frequency fio. Right, data 
taken after calibration; microwaves detuned less than 1 MHz 
from the qubit frequency. 



tally than for peak-finding in high power spectroscopy. 
To increase the |2) state population, the (7r)io pulses (for 
the |0) «-» |1) transitions) are sufficiently short and do not 
use the HD protocol. Using this technique, we measure 
the |1) — )> |2) transition frequency, /21 to within 1 MHz. 



AMPLIFIED PHASE ERROR THEORY 

From numerical simulations, an X^/2 pulse has the 
transformation, (ignoring global phases), given by 
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Using Eqs. (SI) and (S2), we explicitly calculate the 
phase shift for a 2tt rotation coming from four tt/2 pulses. 
For an arbitrary rotation about the x-axis, the gate op- 
eration is 
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FIG. S2. Phase shift versus number n of pseudo- identity op- 
erations. 



Concatenating four positive tt/2 rotations results in 
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where / is the identity. Equation (S6) only acquires a 
global phase. 

Next, we calculate the phase shift for the pseudo- 
identity operation, used in the APE protocol, comprised 
of a positive then a negative = tt/2 rotation. 
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For n applications of the pseudo-identity operation, in 
the limit where 0<e<l, e-^ne 
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and by removing an overall global phase 

^tt/2 ~ (^2e) U = Z2ne • 
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The measured phase shift scales with the number n of 
pseudo-identity operations as shown in Fig. S2. 



Z-PULSE CALIBRATION 

Explicit Z gates are required for the (off-equator) 
Hadamard gate. We calibrate our Z pulse as shown in 
Fig. S3. A static length (full-width at half-max = 6ns) 



with an increasing amplitude Z-pulse is inserted between 
two HD 7r/2 pulses with fixed separation time tfi xe d. The 
probability of measuring the |1) state P\ oscillates with 
increasing Z amp [3]. The arrow indicates the Z amp that 
corresponds to a rotation angle of tt about the z- axis. 
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FIG. S3. Z^ calibration, (a) Left, the Ramsey-type pulse 
sequence to calibrate a Z^ with the X and Y controls using 
the HD protocol described in the manuscript. The sequence 
consists of two, 6ns (FWHM) tt/2- pulses fixed in time with 
tfixed = 24ns and a 6ns (FWHM) Z-pulse centered in between 
them. The separation time is chosen to minimize overlap of 
the pulses. The Z amp increases incrementally. Right, the 
probability of measuring the |1) state Pi as a function of Z- 
pulse amplitude, Z amp . The data are plotted as points with 
best fit as a line. Pi oscillates with increasing magnitude of 
the Z-pulse amplitude. Arrow indicates the Z-pulse amplitude 
equivalent to a tt rotation about the Z-axis. 
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